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i> ■ ABSTRACT 

o 

\ Context. The cut-off frequency of magneto- acoustic gravity (MAG) waves is known to be lowered by the inclined magnetic field. 

This effect allows low-frequency waves to penetrate into the upper atmosphere and the corona. 

Aims. We estimate the modification of the cut-off frequency at various heights and reconstruct the inclination of the magnetic field 
according to the magnetoacoustic-gravity wave theory in a stratified solar atmosphere permeated by a uniform magnetic field. 
Methods. We used the Pixelised Wavelet Filtering (PWF) method to compute narrow-band power maps of SDO/AIA imaging datasets 
in the 1700 A, 1600 A and 304 A bandpasses that correspond to different heights. The cut-off frequency was defined as contours where 
the spectral power droped to the median level. It was measured as a function of the spatial location. We inferred the magnetic field 
inclination according to the MAG wave theory in the low-/? limit and compared it with the potential field extrapolation. 
Results. We analysed intensity oscillations in a symmetric sunspot AR11131 (08 Dec 2010) and an asymmetric sunspot AR11330 
O | (27 Oct 201 1). We reconstructed the magnetic field inclination in the radial direction for the symmetric sunspot and in both radial and 

I ' azimuthal directions for the asymmetric sunspot. 

Conclusions. We observed 3D variation of the main oscillation periods in sunspots. We found that shorter-period oscillations were 
mostly constrained in sunspot umbrae, while longer-period oscillations formed an annular shape enclosing the umbra. Longer periods 
are found to be distributed further away from the sunspot centre. Our results indicate that 3-min oscillation are generated in the 
chromosphere, possibly by the acoustic resonator model, while 5-min and longer-period oscillations seemed to originate in a level 
under the photosphere. The reconstructed field inclinations gives the values of the field inclination that are systematically larger than 
the values obtained by the potential field extrapolation. The inclined magnetic field line can account for 60-80% of cut-off frequency 
lowering only. 
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>' 1. Introduction vo = yg/4nc s = 5.2 mHz (with the corresponding cut-off pe- 22 

riod po = 3.2 min), where g = 274 m/s is the gravitational ac- 23 

~ rA ~^ . t celeration, and c s = 7km/s is the sound speed estimated for 24 

Cj Mangetoacoustic gravity (MAG) waves in the solar atmo- ica] chromospheric conditions . i n the low-/? approximation 25 

ry t sphere are subject to strong dispersion: their properties de- ^ <<; j in a s t 0f coronal actiye ionX ^ cut _ off 2g 

^ pend strongly on their frequency. In particular the parame- f is modified b the etic inclination Vac = VQ cos 27 

•rt ters of the solar atmosphere determine the MAG cut-off fie- (Bel &L 1977) 1 . This helps explaining the existence of low- 28 

quency. MAG waves with frequencies below the cut-off value f waves in the coronaj which ^ bdieved t0 a 29 

M| become evanecent, and therefore cannot reach the upper lay- j Jon of e flux than ^ high . frequency counterpart 30 

* ers of the solar atmosphere The cut-off value also determines ( Fontenla et ^ m3 . Jefferies et aJ 20Q6) A1 ^ stud 31 

9 the oscillating frequency of MAG waves as the response of Qf ^ effect Qf ^ cut _ off fr Qn essive waves ob . 32 

10 the atmosphere to a broadband e.g^ impulsive) excitation (e.g. seryed in ^ solar atm here ides us with a seismo i ogica i 33 

I I Suematsu et al. 1982; Botha et al. 201 1) The cut-off frequency too] for ^ determination of the local t of the etic 34 

1 2 lor MAG waves was derived lor a stratified isothermal solar at- ^ ^ 35 

13 mosphere permeated by a uniform magnetic Md (Bel & Leroy Accordi t0 Bel & L (1977)> the i ow . freque ncy (long- 36 

14 1977; Zhugzhda & Dzhahlov 1984). It generally depends on the riod ^ ^ ^ tQ ^ ^ ^ J atm0 s P herl 37 

15 local plasma fi (ft = l^p/B = 2c JyV A , where p is the gas &[ ^ ^ concentrationS) ots which 38 

16 pressure, llq is the magnetic permeability m vacuum space, B , . , T . , A , u f , , , u 

\ t . o i , 1 1 rr • , i i A ir , , work as MAG waveguides. It was suggested that, due to the 39 

17 is the magnetic field strength, Va is the local Alfven speed, 

18 c s is the local sound speed, and y = 5/3 is the adiabatic in- i The derivations in Bel & Leroy (1977) contain obvious misprints, 

19 dex) and the magnetic field inclination 0. In the high-/? re- a follow-up theoretical study will redo the derivation. However the nu- 

20 gions (J3 » 1, e.g. the photosphere and chromosphere of the merical results and figures are reliable. In this work, we use the trust- 

21 quiet Sun), it is reduced into the pure acoustic cut-off frequency worthy extreme case of p « 1 
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1 effect of the magnetic inclination, the 5-min oscillations ap- 

2 pear in the chromo spheric spicules (De Pontieu et al. 2004) 

3 and in coronal loops at active regions (De Moortel et al. 2002; 

4 De Pontieu et al. 2005). de Wijn et al. (2009) observed the prop- 

5 agation of solar global /?-modes to the chromosphere through the 

6 inclined magnetic fields at the periphery of plage regions. Other 

7 observations of long-period oscillations in the corona could be 

8 attributed to channelling by the inclined field as well (see e.g. 

9 Wang et al. 2009; Marsh et al. 2009; Yuan et al. 201 1). 

10 A direct observation of the magnetoacoustic cut-off in a 

1 1 sunspot has been performed with full Stokes (I, U, V, Q) spec- 

12 tropolarimetry (see Fig. 5 in Bloomfield et al. 2007). The au- 

1 3 thors found that the cut-off frequency closely followed the mag- 

14 netic field modification in the strong-field limit (J3 « 1) of 

15 Bel & Leroy (1977). However no further examination was done 

16 to compare the information with the magnetic inclination in- 

1 7 verted by full Stokes observables. It appeared that the correlation 

1 8 and phase difference of the time series of the intensity variation 

1 9 obtained at various heights of the solar atmosphere looked more 

20 realistic when the offsets at the spatial location were corrected by 

21 considering the magnetic inclination (Bloomfield et al. 2007). 

22 Mcintosh & Jefferies (2006) studied the travel time of narrow- 

23 band signals around a sunspot and found good consistency with 

24 the prediction of Bel & Leroy (1977) for both the quite sun 

25 (J3 » 1) and a sunspot (J3 < 1). Tziotziou et al. (2006) applied the 

26 empirical formula v pe ak(0) ^ l-25v ac (0), where v pea k denotes 

27 the peak frequency at each location, v ac indicates the correspond- 

28 ing magnetoacoustic cut-off frequency, to the waves observed in 

29 a sunspot' s chromosphere, and estimated the magnetic field in- 

30 clination = arccos[v pea k/(1.25 x v )] (Bogdan & Judge 2006). 

31 Reznikova et al. (2012) compared the spectra and phase rela- 

32 tions of the UV and EUV emission intensity at various heights of 

33 the solar atmosphere to identify the features of upwardly propa- 

34 gating waves and found that the variation of the cut-off frequency 

35 across the umbra was consistent with Bel & Leroy (1977). The 

36 follow-up study (Reznikova & Shibasaki 2012) compared the 

37 observational MAG cut-off frequencies in AIA 304 A with the 

38 values obtained by the potential field extrapolation (Sakurai 

39 1982) and found good agreement. 

40 In Bel & Leroy (1977), an adiabatic plasma was implicitly 

41 assumed, while in the photosphere and chromosphere, radiative 

42 loss is of significance. Centeno et al. (2006, 2009) used a lin- 

43 ear wave equation with a radiative cooling term and explained 

44 well the observed phase delay and wave amplitude variation with 

45 height in both sunspots and pores. A detailed study using mul- 

46 tiple spectral lines formed at different heights was presented in 

47 Felipe et al. (2010). The connectivity of different layers of the 

48 atmosphere determined by the phase difference and power am- 

49 plification was demonstrated. Therefore, the study of the cut- 

50 off frequency can also retrieve the plasma parameters associated 

51 with radiative loss. 

52 The cut-off frequency detected in a sunspot atmosphere is a 

53 good indicator of local plasma parameters, therefore it can be 

54 used to diagnose the local atmosphere, e.g. temperature, mag- 

55 netic inclination. In our study, we present the observation of the 

56 spatial distribution and height variation of the cut-off frequency 

57 with SDO/AIA (Lemen et al. 201 1; Boerner et al. 201 1) and in- 

58 f erred the magnetic field inclination. The observational results 

59 were used for the reconstruction of the 3D geometry of the mag- 

60 netic field. The information was deployed to diagnose the mag- 

61 netic field inclination. In this paper, we present the analysed data 

62 sets in Sec. 2. The methods are summarised in Sec. 3. Then we 

63 compare the diagnostic results with the potential field extrapola- 

64 tion in Sec. 4, and conclusions are given in Sec. 5. 



2. Observations 65 

We selected two large sunspots associated with well-developed 66 

active regions in the corona, where outwardly propagating EUV 67 

disturbances were persistently detectable. We chose good obser- 68 

vation intervals without disruptions from either flares or other 69 

transient events. The observation intervals correspond to the 70 

days when the sunspots were crossing the central meridian: they 71 

were least affected by Wilson depression (Loughhead & Bray 72 

1958) and were well exposed for imaging. 73 

The first sunspot was associated with AR1 1131 situated near 74 

the 30° latitude in the northern hemisphere. It crossed the cen- 75 

tral meridian on 08 Dec 2010. It consisted of a strong magnetic 76 

concentration of the south polarity, while the north polarities 77 

were spread sparsely to the west and north of the central sunspot, 78 

therefore the sunspot had a pretty symmetric shape without any 79 

trailing sunspot, see Fig. 1 (top). We used AIA data sets of 1600 80 

A, 1700 A and 304 A from 02:30 to 03:30 UT on 08 Dec 2010. 81 

The cadence time was 24 s for 1600 A and 1700 A images, 82 

and 12 s for 304 A images. This sunspot was also analysed in 83 

Reznikova et al. (2012); Reznikova & Shibasaki (2012). 84 

The second sunspot chosen for this study was the leading 85 

main sunspot on the eastern part of active region AR1 1330 situ- 86 

ated slightly to the north of the solar equator. The trailing sunspot 87 

group of north polarity was about 150 Mm to the main sunspot. It 88 

formed into a non- symmetric shape with a protrusion due to the 89 

tearing of the other polarity (see #los m Fig. 1 (e)). The shape 90 

was inherited at the temperature minimum (1700 A, see Fig. 1 91 

(f) ), the photosphere and transition region (1600 A, see Fig. 1 92 

(g) ). In the chromosphere (304 A, Fig. 1 (h)), the shape was 93 



less pronounced and was replaced by a far-extending fan struc- 94 

ture. The study of propagating EUV disturbances in the coro- 95 

nal fan over this sunspot was presented in Yuan & Nakariakov 96 

(2012). We used one-hour data sets from 04:00 to 05:00 UT on 97 

27 Oct 2011 in 1700 A, 1600 A and 304 A. The instrumenta- 98 

tion and data preparation are presented in Yuan & Nakariakov 99 

(2012). The cadence time was 24 s and 12 s for EU and EUV 100 

images, respectively. 101 

3. Methods 102 

3. 1 . Power map by Pixelised Wavelet Filtering 1 03 



The Pixelised Wavelet Filtering (PWF) method was developed 1 04 

by Sych & Nakariakov (2008). It is based on the wavelet trans- 105 

form and is well validated in determining the spatial, temporal 1 06 

and phase structure of oscillations in an imaging cube I^x, y,t), 1 07 

where A is the wavelength of the data channel, x, y and t are the 1 08 

dicrete spatial locations and the measuring times, respectively. It 1 09 

gives the spatial distribution of the amplitude (power), frequency 1 1 

and phase of a signal in the spectral interval of interest. 1 1 1 

For data sets h(x,y, t), we performed PWF analysis and 112 

obtained narrow-band power maps P* (x, y) (or P A p (x, y) , where 1 1 3 

p = v _1 = 2, • • • , 20min with the resolution of dp = 0.1 min). 114 

In a sunspot, the power distribution normally forms a filled disk 1 1 5 

(short periods / high frequencies) or a power ring (long pe- 116 

riod / low frequencies) concentric at the sunspot centre (e.g. 117 

Sych & Nakariakov 2008; Reznikova et al. 2012). Therefore, to 118 

facilitate quantitative analysis, the power maps were transformed 1 1 9 

into a polar coordinate system with the origin coinciding with 120 

the sunspot centre P* (r, #) , where r is the distance to the sunspot 1 21 

centre, 6 is the polar angle relative to the horizonal line pointing 1 22 

to the solar west (see Fig. 1 (a) and (e)). 123 
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Fig. 1. Multi-instrumental imaging of AR11131 (top, 08 Dec 
2010) and AR11330 (bottom, 27 Oct 2011) and the underly- 
ing sunspots at different heights, (a) and (e): The LOS mag- 
netic field strength determined with HMI. The polar coordinate 
systems used in this study are overlaid, (b) - (d) and (f) - (h): 
AIA intensity images on logarithmic scales at different obser- 
vational wavelengthes show different levels of the sunspots and 
the associated active regions. In order of increasing heights of 
the observed levels, the displayed images illustrate the tempera- 
ture minimum level (1700 A), upper photosphere and transition 
region (1600 A), and the chromosphere (304 A). The dotted lines 
show the border of the sunspot umbra and penumbra determined 
with the 4500 A image. 



1 For the power maps P v (r, 6), the denotation A is dropped to 

2 refer to the general case that is applicable to all wavelengthes. 

3 We computed the maximum power at each location P max (r,6) 

4 (see panel (a) in Fig. 3, Fig. 4 and Fig. 5). The corresponding 

5 peak period is p pQSi k{r,0) (see panel (b) in Fig. 3, Fig. 4 and 

6 Fig. 5). The dependence of /? pea k(^#) shows the spatial distri- 

7 bution of the period with the dominant oscillating power, while 

8 Pmax(r,6) illustrates the significance of these oscillations. The 

9 variance of the spectral power over the spectrum for each pixel 

10 was also calculated, Var(r, 6) = Var(P v (r, 6)) v 2 (see Fig. 3 (d)). 

1 1 This value is a good indicator of the significance of the signal at 

12 a specific location. To estimate the l<x noise level, we took the 

13 average over the spectrum <x(r, 6) = (P v (r, 0)) v (see example in 

14 Fig. 3(e)). For the analysis of the symmetric sunspot, a further 

1 5 averaging over the polar angle 6 was performed to reduce the 6 

16 dimension. Table 1 summarises the denotations and their physi- 

1 7 cal meanings in both ID and 2D cases. 

18 Several typical narrow-band power maps obtained in 304 

19 A bandpass on 08 Dec 2010 are shown in Fig. 2 (also see 

20 Fig. 7 of Reznikova et al. 2012). For the periods less than 

21 3.2 min (frequencies geater than 5.2 mHz), the dominant power 

22 fills up the umbral region. The inhomogeneity of the power dis- 

23 tribution implies a fine structure in the umbra (see the study 

24 in Jess et al. 2012, using imaging data with a resolution of 

25 50 km/pixel in contrast to 430 km/pixel of AIA images). For 

26 the periods greater than 3.2 min, the power of significant os- 

27 cillations is normally concentrated in an annular structure en- 

28 closing the umbra. The annulus expands with the increasing 



2 The subscript v means the operation(s) along the spectral dimen- 
sion. With the same philosophy, the subscripts 6 and r mean the opera- 
tion^) along the polar angle dimension and radial distance dimension, 
respectively. 



period. We interpret this phenomenon as modification of the 29 

cut-off frequency by the inclined magnetic field (Bel & Leroy 30 

1977; De Pontieu et al. 2004, 2005; Mcintosh & Jefferies 2006, 31 

see discussions in Reznikova et al. (2012)). As the magnetic 32 

field lines approach the outer penumbra, the inclination angle 33 

becomes larger, therefore the cut-off frequency is lowered to a 34 

smaller value (Bel & Leroy 1977), and allows the channelling of 35 

upwardly propagating waves of lower frequency (longer period) . 36 

3.2. 1D reconstruction in symmetric sunspot (08 Dec 2010) 37 

3.2.1 . Reconstruction of field line inclination by MAG cut-off 38 

For the symmetric sunspot observed AR11131 (08 Dec 2010), 39 

we present the example of 304 A data in Fig. 3. The distribution 40 

of the peak power P max (r, 6) and the corresponding peak period 41 

/Vak(r, 6) are illustrated in Fig. 3 (a) and Fig. 3 (b), respectively. 42 

The dominant oscillating power fills in the umbra and gradually 43 

decreases further out in the penumbra and arcade region. In the 44 

umbra, the peak period is the well-known chromospheric 3 -min 45 

oscillation (see review paper Bogdan & Judge 2006). It fills up 46 

the whole umbra and a tiny part of the penumbra. In the penum- 47 

bra, the peak period gradually increases with the radial distance 48 

from sunspot centre, indicating the modification of the cut-off 49 

frequency by inclined magnetic field lines. The variation with 50 

the polar angle is found to be very small. We averaged the power 51 

along the polar angles and obtained the normalised power P v (r) 52 

as a function of the distance to the sunspot centre, r, and the fre- 53 

quency v, as shown in Fig. 3 (c). We defined the contour of the 54 

global median in P v (r) as the curve of the cut-off period p ac (r) 55 

(or frequency v ac (r) = l/p ac (r), see Fig. 3 (c)). 56 

According to Bel & Leroy (1977), the MAG cut-off fre- 
quency in the low-/? approximation is modified when the mag- 
netic field line is inclined from the vertical direction, v ac = 
vo cos cp. We estimated the inclination angle by 

<p(r) = arccos = arccos — — (1) 



vo Pac(r) 

where vo = 5.2 mHz (po = 3.2 min) is the acoustic cut-off fre 
quency (period) (see Fig. 3 (f)). 

3.2.2. Error and significance analysis 

According to (1), the uncertainty of (p is calculated as 

Avac(r) cos 2 Apzc(r) 



A0(r) 



yl - Vac 2 0) 



y/\ - COS 2 P0 



(2) 



which requires the estimation of A/? ac . This quantity can be as- 
sumed to be proportional to the noise level <x(r), and the spread 
(standard deviation) of the spectral power over the polar angle 
s(r). Ap ac should also be inversely proportional to the maximum 
spectral power P max (r) and the square root of the number of the 
data point in polar angles, y/Ne. The noise level cr(r) is found to 
vary between 0.16 and 0.36. When we normalise it with its total 
median, it gives a factor between 0.52 and 1.13 to the uncer- 
tainty. We did the same to the spread s(r) and obtained another 
factor between 0.31 and 1.4 added to the error bar. The value of 
Pmax(r) changes from 1 to about 0.4 with r, giving a factor from 
about 1 to about 2.5 to the total error bar. We choose a constant 
dp = 0.1 min, the spectral resolution of PWF analysis, therefore 
the total uncertainty becomes 

a n = °"( r ) ^( r ) d P n x 

Pac[n MedHr)) Med(*(r)) y/N e P^(r) ' K } 
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Table 1. Notations used in the paper. 
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Physical meaning 


2D denotation 


Physical meaning 


ID denotation 


2D spectral power map 


P v (r,6) 




ID spectral power map 


P v (r) = (P v (r,6)) d 


2D spectral noise level 


<r(r,Q) 


= <Pv(r,0)\ 


ID spectral noise level 


<r(r) = <P v (r)> v 


2D spectral power variance 


Var(r,<9 


) = Var(P v (r,(9)) v 


ID spectral power variance 


Var(r) = Var(P v (r)) v 


2D spectral peak power 


Pmdodj" •> 


0) = max(P y (r,(9)) y 


ID spectral peak power 


^maxO) = max(P v (r)) v 


2D peak period 




0) 


ID peak period 


PpeakM for all 


2D MAG cut-off period 


p ac (r,0 


) 


ID MAG cut-off period 


p ac (r) for all 6 


2D MAG cut-off frequency 


Vac (r, 6) 


= Vp ac (r,6) 


ID MAG cut-off frequency 


VacO) = l/p ac (r) 








ID power spread 


s(r) = stdev(P y (r,6>)) | v=Vac 
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Fig. 2. Typical narrow band power maps of 304 A bandpass in polar coordinate made for 08 Dec 2010. 



1 When the cut-off period approaches the detection limit (20 

2 min in our application), the observation of the cut-off effect be- 

3 comes less reliable. Applying this limit, we estimated the re- 

4 gions of significant measurements. We estimated the noise in 

5 the spectra as cr(r,6) = (P v (r,6)) v . The noise level is illus- 

6 trated in Fig. 3 (e). It was further averaged over the polar angles, 

7 <x(r) = (P v (r, 6)) e y . The measurements with the noise level dis- 

8 tributed above 68% (more than a standard deviation above the 

9 mean) of the noise are considered as less reliable (see the cross- 

10 hatched regions in all panels of Fig. 3). An alternative way is 

1 1 to consider the variance of the power maps, which reflects the 

12 contrast of the oscillating powers, giving us the significance of 

13 the measurements. We excluded 23% (lower than one standard 

1 4 deviation below the mean) of the smaller portion of the power 

15 variance distribution. As the results determined with the use of 

16 either power variance or noise distribution are almost similar, 

1 7 we utilised only the results quantified by the noise level in the 

18 significance analysis. 



19 3.3. 2D reconstruction in asymmetric sunspot (27 Oct 2011) 



20 In 2D spectral power maps, we consider the areas with the spec- 

21 tral power above the median of the corresponding narrow-band 

22 power map as regions that allowed upward propagating MAG 

23 waves. The borders of these areas were defined as contours of 

24 the MAG cut-off. By labelling the cut-off contours, we were able 

25 to reconstruct the field line inclination. 



The reconstruction steps are similar to those in the ID case. 
The estimation of the error, A/? ac (r, 6), is simplified, 

cr(r, 0) 



A/?acO, 



dp 



(4) 



Med(<r(r,0)) P v (r,6) 

We masked out the regions with the less significant measure- 26 

ments. The masks were designed with two criteria: 1) the pixel 27 

should be within the penumbra and umbra; 2) the noise level 28 

in that pixel should not exceed the 68% level in <x(r, 6) distri- 29 

bution. Only the measurements within the mask are shown in 30 

Fig. 6, Fig. 7 and Fig. 8 31 



3.4. Magnetic field extrapolation 



32 



We compare the field reconstruction with 3D magnetic field ex- 33 

trapolation. Potential (current- free) fields were extrapolated us- 34 

ing the MAGPACK2 package (Sakurai 1982). The current-free 35 

model in this code is based on the Green ' s function method and 36 

requires a line-of- sight magnetogram as the bottom boundary 37 

condition. Although, we can not a priori exclude the presence 38 

of the electric currents in the considered volume, we have cho- 39 

sen this model because of the following reasons: 1) both of ARs 40 

did not produce any flares within or close to our observational 41 

intervals; 2) we are interested in field structures in the layers be- 42 

low the corona where the dense material can influence the field 43 

inclination. Therefore, the choice of the non-potential field is un- 44 

likely to produce better approximation at these heights. 45 

Fig. 9 presents a 2D projection of reconstructed magnetic 46 

field lines overlaid on the photospheric line-of-sight magne- 47 

togram (a) sunspot AR11131 (08 Dec 2010) and (b) sunspot 48 



4 



Ding Yuan et al.: MAG cut-off 




Polar angle (71) Polar angle (71) Inclination (degree) 

Fig. 3. The oscillating power information extracted from the 304 A data set of the sunspot AR11131 (08 Dec 2010), shown in the 
polar coordinates (r and 6). The dashed lines mark the borders of the umbra and penumbra. The unreliable region is cross-hatched, 
a) The peak power distribution, b) The peak period distribution, c) The ID spectral power map as a function of r and p. The 
contour in asterisk shows the cut-off period, d) The spectral power variance distribution, e) The spectral noise distribution, f) The 
reconstruction of the magnetic field inclination. 



1 AR1 1330 (27 Oct 201 1). The colour coding represents the mag- 

2 netic field strength, using blue for southern (negative) and red 

3 for northern (positive) polarities. For the sunspot AR1 1 131, we 

4 selected sub-area of 160" x 160" from the SOHO/MDI full-disk 

5 magnetogram taken at 04:51 UT, as shown in Fig. 9(a). For the 

6 sunspot AR1 1330, the HMI line-of-sight magnetogram was ob- 

7 tained as a level- 1.5 FITS file. A larger size of 300" x 300" was 

8 selected to include the sufficient amount of positive polarity fol- 

9 lowing the main sunspot. The magnetic field was computed with 

10 a 2" mesh spacing using a spherical boundary. 

11 4. Results 

12 4.1. Spatial and height distribution of dominant oscillations 

13 Based on the ID power distribution, panel (c) of Fig. 3, Fig. 4 

14 and Fig. 5, it is clearly seen that 3-min oscillations are very 

15 strictly constrained within the umbra region, while 5-min oscil- 

16 lations dominate only at a thin annulus (about 2.9 - 5.0 Mm in 

1 7 width) between the umbra and penumbra. Although 3-min oscil- 

1 8 lation is strongest in the umbra, it is predominant within the disks 

19 of radii about 2.1, 5.8 and 7.2 Mm in the 1700 A (temperature 

20 minimum), 1600 A (upper photosphere) and 304 A (chromo- 

21 sphere) data, respectively. Although 3-min oscillation can prop- 

22 agate through both vertical and inclined magnetic field, it is con- 

23 strained only in the umbra. This implies that the 3-min oscilla- 

24 tion is not likely to originate from the solar interior. Instead, it 

25 is consistent with the acoustic resonator model (Zhugzhda 2008; 

26 Botha et al. 201 1). The 5-min oscillation is well pronounced in 

27 the 1700 A data. It appears that there is strong interaction be- 



tween the magnetic fields and the ambient plasma enclosing the 28 

sunspot umbra, where the magnetic field is medium in strength 29 

and slightly inclined and has a strong longitudinal gradient. It 30 

seems there is strong p-mode absorption at the umbra border. 31 

The maximum energy transmission was predicted to occur at the 32 

height where the Alfven speed is equivalent to the local sound 33 

speed and at an attack angle of about 30° (Cally et al. 2003; 34 

Schunker & Cally 2006). Thus the 5-min oscillation may form 35 

at the temperature minimum (about 500 km above the photo- 36 

sphere). The 2D power maps also confirm these results, see pan- 37 

els a) and b) in Fig. 3, Fig. 4, Fig. 5, Fig. 6, Fig. 7 and Fig. 8. 38 

In panel (a) of Fig. 3 and Fig. 6, the power of 3-min oscillations 39 

show variations within the umbra. This implies fine structuring 40 

inside the umbra. Stronger sources of oscillations were found 41 

better correlated with the 3-min propagating disturbances in the 42 

corona loops in 171 A bandpass (Jess et al. 2012; Sychetal. 43 

2012). A dedicated phase study of the oscillations at various 44 

heights of the sunspot atmosphere will be presented in a follow- 45 

up paper. 46 

4.2. Correlation of the peak and cut-off frequencies 47 

The detected cut-off frequencies are plotted versus the corre- 48 

sponding peak frequencies, see Fig. 10. The correlation between 49 

them is very high: the Pearson's correlation coefficients R are 50 

0.98, 0.94 and 0.99 for the 1700 A, 1600 A and 304 A data, 51 

respectively. We fitted a linear relationship (v ac = av pea k + b) 52 

to the data using the MPFITEXY routine 3 (Williams et al. 53 

3 http://purl.org/mike/mpfitexy 
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Fig. 4. The same as Fig. 3 but in 1600 A 



1 2010). The MPFITEXY routine depends on the MPFIT pack- 

2 age 4 (Markwardt 2009). All the three lines are very close to one 

3 another, the slopes are 0.80 ± 0.01 5 , 0.75 + 0.01 and 0.74 + 0.01 

4 for the 1700 A, 1600 A and 304 A data, respectively. The re- 

5 suits agree with the empirical relation (v pea k — 1.25v ac ) between 

6 the peak and cut-off frequency as suggested in Bogdan & Judge 

7 (2006); Tziotziou et al. (2006), but we obtained non-zero inter- 

8 cepts in the linear fits, -(0.78 + 0.01) mHz , -(0.47 + 0.01) mHz 

9 and -(0.78 + 0.01) mHz for 1700 A, 1600A and 304 A, respec- 

1 tively (see Fig. 10). Large deviations from the linear model occur 

1 1 at Vpeak — 3 mHz (5-min band), strongest in 1700 A data, and also 

1 2 visible in 1600 A and 304 A as well. 

13 4.3. Reconstruction of the magnetic field inclination 

14 As described in Sec. 3.2, the results of ID reconstruction are il- 

1 5 lustrated in Fig. 3 (304 A), Fig. 4 (1600 A) and Fig. 5 (1700 A). 

1 6 Comparison of our results obtained with the seismological tech- 

17 nique with the potential field extrapolation is shown in Fig. 11. 

1 8 The inclination angles were retrieved at 500 km (1700 A), 1 100 

19 km (1600 A) and 2200 km (304 A) above the solar surface. 

20 The magnetic inclination was estimated with both the theoretical 

21 acoustic cut-off frequency vo = 5.2 mHz and the corresponding 

22 observationally determined cut-off frequency v{J, where the in- 

23 dex A corresponds to the observational channel used in its deter- 

24 mination. 

25 Fig. 1 1 shows clearly that the general profiles of field incli- 

26 nations reconstructed with the MAG cut-off at 1600 A, 1700 A 

27 and 304 A, agree very well with the potential field extrapolation 

28 in the region from 0.2r p to 0.8r p , where r p is the average radius 

4 http://cow.physics.wisc.edu/ craigm/idl/idl.html 

5 The fit errors less than 0.01 were considered as 0.01. 



of the penumbra. The average offsets are 32.7/28.8, 27.9/26.6, 29 

30.9/29.7 for v /v* reconstructions in 1700 A, 1600 A and 304 30 

A, respectively. The results are slightly improved by using the 31 

observed maximum cut-off frequency, but the large offsets can- 32 

not be fully removed. This implies that the inclined magnetic 33 

field is not the only factor that influence the MAG cut-off fre- 34 

quency, it only accounts for 60-80% of the lowering in the cut- 35 

off frequency. In the inner umbra (from to 0.2r p ), the devia- 36 

tions from the potential field extrapolation are about 20° — 30°. 37 

These deviations may arise from using a constant vo; choosing 38 

the observed vjj for the corresponding bandpass appears more 39 

realistic. We obtained a flat distribution in all bandpasses in the 40 

inner umbra. It implies that the wave characteristics inside the 41 

inner umbra may be different from other regions of the sunspot, 42 

or that the impact of the plasma temperature needs to be con- 43 

sidered as well. We notice that at 0.8r p -1.2r p in the 1600 A and 44 

1700 A data, a region of low inclination was obtained. In our 45 

method, cos(+0) = |cos(180° + 0)| can not be differentiated 46 

(the 180° ambiguity problem, see Metcalf et al. 2006). By con- 47 

sidering this effect, this low-inclination region in the penumbra 48 

border can be connected with the return-flux region with the 49 

magnetic vector pointing downwards, instead of upwards. The 50 

green asterisk denotes the return flux in the potential magnetic 51 

field corrected into [0°, 90°]. The results are consistent with the 52 

values estimated with the use of the MAG cut-off, both in the in- 53 

clination and the distance to the sunspot centre. Such behaviour 54 

is consistent with the return-flux sunspot model (Fla et al. 1982; 55 

Osherovich 1982). The return-flux region is found to spread at 56 

about 0.8r p - 1.2 r p . If we assume a parabolic shape of the mag- 57 

netic field lines, the returning flux at about 0.4r p -0.6 r p extends 58 

to the height of about 0.14 r p - 0.26 r p and returns to the surface 59 

about 0.4 - 0.6 r p apart from the source. This effect is less pro- 60 

nounced in the 304 A data, in which the magnetic field shows a 61 
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Fig. 5. The same as Fig. 3 but in 1700 A 



1 steady extension outwards the sunspot. The agreement between 

2 the inferred magnetic inclination and the values obtained with 

3 the potential field extrapolation is better in the return-flux region 

4 compared with the inner sunspot region. The return fluxes were 

5 normally found at the periphery of the sunspot penumbra with a 

6 weaker magnetic field and higher temperature, so we infer that 

7 other effects, e.g. radiative cooling, should be accounted as well 

8 to understand the lowering of the MAG cut-off frequency. 

9 In 2D reconstruction, we compare the reconstructed mag- 

1 netic field inclination, panel (e) of Fig. 6, Fig. 7 and Fig. 8, with 

1 1 the potential field extrapolation shown in Fig. 12. Our method re- 

1 2 constructed the magnetic inclination distribution: both the global 

13 profiles, fine structure inside the sunspot umbra and the protru- 

1 4 sion at the penumbra. In the 1600 A and 1700 A data, the return 

15 flux is reconstructed very well. It appears that in 2D geometry, 

1 6 no apparent physics modified the result, except the morphology. 

1 7 The potential field extrapolation reflects the general profiles, in- 

1 8 eluding the protrusion structure and the return flux. The average 

19 offsets are larger, 41.0, 37.8 and 37.0, which may arise from the 

20 error accumulation, since the outliers can not be effectively ex- 

21 eluded. 



22 5. Conclusions 

23 We observed the 3D period variations of compressive oscilla- 

24 tions in two sunspots. Our results confirm the previous findings 

25 that short-period (high-frequency) oscillations are constrained 

26 inside the sunspot umbra, while long-period (low-frequency) 

27 oscillations are mainly present outside the umbra. Longer the 

28 period is, further away the oscillations are observed from the 

29 sunspot centre. The predominant oscillation inside the umbra is 

30 the well-known 3-min oscillation. The variation over different 

31 heights indicates that the 3-min oscillation is not sourced from 



the solar interior, rather it is more consistent with the acous- 32 

tic resonator model (Zhugzhda 2008; Botha et al. 201 1). The 5- 33 

min oscillations form a ring-shape at the umbra border. We may 34 

speculate that the strong power of 5 -min oscillations may indi- 35 

cate strong interaction of acoustic waves with the magnetic field, 36 

where the local condition favoured the absorption of global solar 37 

/?-modes (Cally et al. 2003; Schunker & Cally 2006). 38 

We seismologically reconstructed the magnetic field inch- 39 

nation with the use of the Bel & Leroy (1977) formalism, and 40 

found its results to be consistent with the results obtained 41 

with potential field extrapolation. It confirmed the conclusion 42 

in Bel & Leroy (1977) that the inclined magnetic field lowers 43 

the cut-off frequency. The MAG-reconstructed magnetic inch- 44 

nations are generally larger than the potential field extrapolation. 45 

The small discrepancy may be attributed to the other effects,e.g. 46 

radiative cooling (Centeno et al. 2006, 2009; Felipe et al. 2010). 47 

The effect of the magnetic field inclination accounts for about 48 

60-80% of the lowering in the cut-off frequency. The seismo- 49 

logical method applied in our study provides an alternative way 50 

to reconstruct the magnetic field. Moreover, it may provide the 51 

boundary or initial conditions for the magnetic field extrapola- 52 

tion. 53 

In our method, we assumed J3 « 1 , which may be incorrect 54 

in some parts of the sunsports. The method can be improved by 55 

considering an observational or empirical ft profiles at various 56 

heights of the sunspots. If the 3D magnetic field vector field is 57 

available, the magnetic field strength and inclination can be cal- 58 

culated, we are able to reconstruct the sound speed and local 59 

temperature profiles in both ID and 2D at various heights. 60 
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Fig. 9. Results of the potential field extrapolation using the MDI magnetogram for sunspot AR1 1 131 on 08 Dec 2010 (a) and with 
the HMI magnetogram for sunspot AR1 1330 on 27 Oct 201 1 (b). 



1 Federation. The AIA data were used as courtesy of the SDO AIA science team. 

2 This research has made use of NASA's Astrophysics Data System. 
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Fig. 6. a) - d) are the same as a), b), d) and e) in Fig. 3. e) and f) Fig. 7. The same as in Fig. 6 but in 1600 A 
are the inclination angle and its uncertainty reconstructed in 2D. 
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Fig. 8. The same as in Fig. 6 but in 1700 A 
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Fig. 10. Correlation between the cut-off frequency and peak fre- 
quency for 1700 A (a), 1600 A (b) and 304 A (c) data. The dia- 
mond symbols represent the measurement at sunspot AR11131 
observed on 08 Dec 2010, while dot symbols indicate the data 
of sunspot AR1 1330 observed on 27 Oct 201 1. The data are fit- 
ted with a linear relationship v ac = av pea k + b, as plotted in red 
solid line. The fit parameters, as well as the Pearson's correlation 
coefficient, are displayed in each panel. 
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Fig. 11. Comparison of magnetic field inclination reconstructed 
by magneto-acoustic cut-off with potential field extrapolation for 
sunspot AR11131 (08 Dec 2010) in 1700 A (a), 1600 A (b) and 
340 A (c) bandpasses. The inclination angles obtained with the 
potential field extrapolation are represented in grey asterisk as 
the scatter plot, while the data in green asterisks denote the ab- 
solute values of inclinations of the return flux. The inclination 
angles reconstructed with vo = 5.2 mHz are plotted in the red di- 
amonds, while those reconstructed with the observational max- 
imum cut-off frequency vj* are plotted in blue diamonds. The 
dashed lines show the border of the sunspot umbra and penum- 
bra, The cross-hatched region marks the measurements below 
the significant level. 
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Fig. 12. Left column: magnetic field inclination obtained with 
vo reconstruction. Right column: the inclination angle obtained 
from potential field extrapolation for AR1 1330 on 27 Oct 201 1. 
From top to bottom rows, the panels correspond to 304 A, 1600 
A and 1700 A data. The dotted lines mark the umbra and penum- 
bra borders. 
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